Detection of Supernova Neutrinos by Neutrino-Proton Elastic Scattering 
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We propose that neutrino-proton elastic scattering, u + p —* v + p, can be used for the detection 
of supernova neutrinos in scintillator detectors. Though the proton recoil kinetic energy spectrum 
is soft, with T p ~ 2E„/M P , and the scintillation light output from slow, heavily ionizing protons is 
quenched, the yield above a realistic threshold is nearly as large as that from u e + p —* e + + n. In 
addition, the measured proton spectrum is related to the incident neutrino spectrum, which solves 
a long-standing problem of how to separately measure the total energy and temperature of v^, u T , 
Vp, and v T . The ability to detect this signal would give detectors like KamLAND and Borexino a 
crucial and unique role in the quest to detect supernova neutrinos. 

PACS numbers: 97.60.Bw, 13.15.+g FERMILAB-Pub-02/087-A 



I. INTRODUCTION 

When the next Galactic supernova occurs, approxi- 
mately 10 4 detected neutrino events are expected among 
the several detectors around the world. It is widely be- 
lieved that these 10 4 events will provide important clues 
to the astrophysics of the supernova as well as the prop- 
erties of the neutrinos themselves. Interestingly, recent 
breakthroughs in understanding solar and atmospheric 
neutrinos each occurred when the accumulated samples 
of detected events first exceeded 10 4 . 

But will we have enough information to study the su- 
pernova neutrino signal in detail? Almost all of the de- 
tected events will be charged-current u e + p — ► e+ + n, 
which will be well-measured, both because of the large 
yield and because the measured positron spectrum is 
closely related to the neutrino spectrum. Because of 
the charged-lepton thresholds, the flavors v^, v T , v^, and 
v T can only be detected in neutral-current reactions, of 
which the total yield is expected to be approximately 
10 3 events. However, as will be discussed below, in gen- 
eral one cannot measure the neutrino energy in neutral- 
current reactions. This paper presents an exception. 
These four flavors are expected to carry away about 2/3 
of the supernova binding energy, and are expected to have 
a higher temperature than v e or v e . However, there is no 
experimental basis for these statements, and at present, 
numerical models of supernovae cannot definitively ad- 
dress these issues either. If there is no spectral signature 
for the neutral-current detection reactions, then neither 
the total energy carried by these flavors nor their tem- 
perature can be separately determined from the detected 
number of events. 

But it is crucial that these quantities be measured. 
Both are needed for comparison to numerical supernova 
models. The total energy is needed to determine the 
mass of the neutron star, and the temperature is needed 
for studies of neutrino oscillations. At present, such stud- 
ies would suffer from the need to make model-dependent 
assumptions. This problem has long been known, but 



perhaps not widely enough appreciated. In this paper, 
we clarify this problem, and provide a realistic solution 
that can be implemented in two detectors, KamLAND 
(already operating) and Borexino (to be operating soon) . 
The solution is based on neutrino-proton elastic scatter- 
ing, which has been observed at accelerators at GeV en- 
ergies, but has never before been shown to be a realistic 
detection channel for low-energy neutrinos. Some of our 
preliminary results have been reported at conferences [|lj . 



II. CROSS SECTION 

The cross section for neutrino-proton elastic scattering 
is an important prediction j!| of the Standard Model, and 
it has been confirmed by extensive measurements at GeV 
energies (see, e.g., Ref. ||). At the energies considered 
here, the full cross section formula ||, || can be greatly 
simplified. At low energies, the differential cross section 
as a function of neutrino energy E u and struck proton 
recoil kinetic energy T p (and mass M p ) is 
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The neutral-current coupling constants between the ex- 
changed Z° and the proton are 
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where the factor 1.27 is determined by neutron beta de- 
cay and its difference from unity is a conseq uen ce of the 
partially conserved axial current. Equation (2.1) may be 
obtained directly by summing the contributions from the 
valence quarks. The cross section for antineutrinos is ob- 
tained by the substitution c A — > —c A . At high energies, 
the primitive couplings are functions of q 2 /M 2 , where 



M ~ 1 GeV (the proton mass or the dipole form-factor 
masses); since q 2 — 2M p T p ~ E 2 , this variation may 
be safely neglected here. At order E y /M p: there is also 
a weak magnetism term which we have neglected. This 
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would appear inside the square brackets in Eq. (2.1) as 



4T p E v c m ca ■ 



(2.4) 



where cm ~ 1.4 depends on the proton and neutron mag- 
netic moments [|] . This term is thus positive for neutri- 
nos and negative for antineutrinos. Besides being numer- 
ically small (less than a 10% correction), this term will 
cancel in the measured differential cross section due to 
the indistinguishable contributions of neutrinos and an- 
tineutrinos. For z/ju and j/ r , we assume the same fluxes 
and spectra for particles and antiparticles (as well as each 
other) ; the weak magnetism term above causes small cor- 
rections to the emitted spectra that we can neglect 
here. For v e and v e , the expected fluxes and spectra are 
different from each other, but at the lower energies of 
these flavors the whole correction is very small. Other 
than the above points, Eq. (2.1) is correct to all orders 
in E v /Mp. As will be emphasized below, our results are 
totally independent of oscillations among active flavors, 
as this is a neutral-current reaction. 

We use the struck proton kinetic energy in the labo- 
ratory frame as our kinematic variable. For a neutrino 
energy E v , T p ranges between and T™ aa: , where 
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The maximum is obtained when the neutrino recoils 
backwards with its original momentum E v , and thus the 
proton goes forward with momentum 2E U . The other 
kinematic variables can be related to T B , and are 



cos 9 p = 
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where 8 p and Q. 



are the angles of the final proton and 
neutrino with respect to the direction of the incident neu- 
trino. In a scintillator-based detector, the proton direc- 
tion cannot be measured, so these expressions are useful 
just for checking the cross section and kinematics. 

If we take (E u —T p ) 2 ~ E 2 (i.e., keeping only the lowest 
order in E u /M p , a very good approximation), then the 
differential cross section is very simple: 
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(2.8) 

Since ca 3> cy, this form makes it clear that the largest 
proton recoils are favored, which is optimal for detection. 
Plots of da/dTp for fixed E v are shown in Fig. [j]. Note 
that these slope in the opposite sense of the corresponding 
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FIG. 1: The differential cross section as a function of T p 
for fixed E v . Note the rise at large T p , indicating that large 
kinetic energies are preferred. From left to right, the lines are 
for E v = 20, 30, 40, 50, and 60 MeV. 



da/dT e curves for — e~ scattering. The difference is 
simply due to the very different kinematics. For neutrino- 
proton elastic scattering, T™ aa: ~ 2E 2 /M p <C E v , while 
for neutrino-electron scattering, T™ ax ~ E v . In this 
limit, the neutrino (ca) and antineutrino {—ca) cross sec- 
tions are identical. If cy is neglected and the differential 
cross section is expressed in terms of cos#„, it follows 
the form 1 — 1/3 cos 9 V expected for a non-relativistic ax- 
ial coupling (i.e., a Gamow- Teller matrix element). The 
total cross section is 



G%E: 



v - (4 + 34) 



(2.9) 



As expected, this is of the same form as the total cross 
section for the charged-current reaction D e + p — > e + + n 
(see, e.g., Rcf. ||). In the neutral-current case, the vector 
coupling nearly vanishes, and the axial coupling is half 
as large as in the charged-current channel, making the 
total cross section approximately 4 times smaller. This 
factor of 4 can be immediately obtained by considering 
the product of the couplings and the propagator factor, 
and using the definition of 9w- 

It is also interesting to compare the neutrino-proton 
elastic scattering cross section with that for neutrino- 
electron elastic scattering (for so that only the neutral- 
current part is compared). Again, the different kinemat- 
ics, reflected in the maximum kinetic energies, are cru- 
cial. The cross section for neutrino-electron scattering is 
much smaller: 



<rtot(vn 



vtotiv^ +p) 



G 2 F E 2 



(2.10) 



which is ~ 10~ 2 for our range of energies 1 . 

In the above expressions, we have neglected contribu- 
tions from strange sea quar ks in the proton |Q . Strange- 
quark effects can enter Eq. (|2.l|) in three ways M. First, 
the vector form factor cy is modified by the strangeness 
charge radius squared (r s ) 2 by a term proportional to 
q 2 (r s ) 2 . Since our q 2 is so low, this is negligible. Second, 
the magnetic form factor cm is modified by the strange 
magnetic moment of the proton /j, s . This is numerically 
small, and appears only in the small weak magnetism 
correction (see above). Third, the strange-quark contri- 
bution As to the nucleon spin gives an isoscalar contri- 
bution to the axial form factor ca, as 



ca 



ca 



1.27 



As 



(2.11) 



The rule for the cross section given above, of using ca 
for neutrinos and —ca for antineutrinos, is also true for 
the combined ca expression given here || . The value of 
As is very poorly known from experiment, and is perhaps 
As = — 0.15±0.15 [||. Since ca S> cy, this could increase 
the differential cross section by approximately 30%, with 
an uncertainty of the same size. It is important to note 
that the As contribution would not change the shape of 
the differential cross section, since ca 3> cy . It may 
be possible to measure As directly via neutrino-proton 
elastic scattering at ~ 1 GeV in MiniBooNE ||. 



III. SUPERNOVA NEUTRINOS 

In this paper, we characterize the supernova neutrino 
signal in a very simple way, though consistently with nu- 
merical supernova models |jl0|| . The change in gravita- 
tional binding energy between the initial stellar core and 
the final proto-neutron star is about 3 x 10 53 ergs, about 
99% of which is carried off by all flavors of neutrinos and 
antineutrinos over about 10 s. The emission time is much 
longer than the light-crossing time of the proto-neutron 
star because the neutrinos are trapped and must diffuse 
out, eventually escaping with approximately Fermi-Dirac 
spectra characteristic of the surface of last scattering. In 
the usual model, v^, v T and their antiparticles are emit- 
ted with temperature T ~ 8 MeV, v e has T ~ 5 MeV, and 
v e has T ~ 3.5 MeV. The temperatures differ from each 
other because v e and v e have charged-current opacities 
(in addition to the neutral-current opacities common to 
all flavors), and because the proto-neutron star has more 
neutrons than protons. It is generally assumed that each 
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FIG. 2: The relative spectra of neutrinos that interact via 
neutrino-proton elastic scattering. From left to right in peak 
position, the curves correspond to u e , v e , and the sum of z/ M , 
v T , and v T . The flux factors N v = {E B /G)/{E V ) ~ l/T 
have been included in the weighting. 



of the six types of neutrino and antineutrino carries away 
about 1 /6 of the total binding energy, though this has an 
uncertainty of at least 50% |ll| . The supernova rate in 
our Galaxy is estimated to be (3 ± 1) per century (this is 
reviewed in Ref . j^] ) . 

The expected number of events (assuming a hydrogen 
to carbon ratio of 2 : 1) is 
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(3.1) 



(Though written slightly differently, this is equivalent to 
the similar expression in Ref. p3|.) We assume D = 10 
kpc, and a detector fiducial mass of 1 kton for Kam- 



LAND. As written, Eq. (3.1) is for the yield per flavor 



assuming that each carries away a portion E of the to- 
tal binding energy (nominally, Eb — 3 x 10 53 ergs, and 
E = Eg/6). The thermally-averaged cross section (the 
integral of the cross section with normalized Fermi-Dirac 
distribution) is defined for each CH2 "molecule", and a 
factor of 2 must be included for electron or free proton 
targets. The spectrum shape of supernova events which 
interact in the detector is given by the product of the 
cross section and a Fermi-Dirac distribution, i.e., 



It is interesting to note that the total rates of neutrino-proton 
elastic scattering events from solar neutrinos are huge: in the 1 
kton KamLAND detector, the rates from the pp, 7 Be, and 8 B 
fluxes are very roughly 10 3 /day, 10 3 /day, and 10 2 /day, respec- 
tively; however, these are only at very low (unquenched) proton 
kinetic energies of approximately 0.2 keV, 2 keV, and 200 keV, 
respectively. 



(F \ E » 
dE v ~ a[ v) l + exp(E v /T) 



(3.2) 



For a cross section a ~ E 2 , this peaks at about 4T (for 
comparison, the average neutrino energy before weighting 
by the cross section is 3.15T), and the yield N ~ T. 
Prior to this paper, the largest expected yield in any 



oil or water detector was from v e +p — > 



As noted 
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FIG. 3: The thermally- averaged differential cross section for 
Fermi-Dirac distributions of temperature T — 3.5, 5, 8 MeV, 
from left to right. This illustrates how the proton spectrum 
changes with the assumed neutrino temperature (since this is 
a neutral-current cross section, it is flavor-independent). 
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FIG. 4: The true proton spectrum in KamLAND, for a stan- 
dard supernova at 10 kpc. In order of increasing maximum 
kinetic energy, the contributions from u e , v e , and the sum of 
v^, v T , P m , and v T are shown with dashed lines. The solid 
line is the sum spectrum for all flavors. Taking the detector 
properties into account substantially modifies these results, 
as shown below. 



in Section II, the total cross sections for charged-current 
v e + p — > e + + n and neutral-current v + p — > v + p 
have similar forms, though the latter is about 4 times 
smaller. However, this is compensated in the yield by 
the contributions of all six flavors, as well as the higher 
temperature assumed for and v T (T = 8 MeV instead 
of 5 MeV). Thus, the total yield from v + p ^ v + p is 
larger than that from v e + p — ► e + + n, when the detector 
threshold is neglected. 

Taking into account radiative, recoil, and weak mag- 
netism corrections, the thermally-averaged cross section 
for v e +p -> e+ + n at T = 5 MeV is 44 x 1CT 42 cm 2 (for 
2 protons) These corrections reduce the thermally- 
averaged cross section by about 20%, and also correct the 
relation E e — E v — 1.3 MeV. The total expected yield 
from this reaction is thus about 310 events in 1 kton. 

In Fig. |2|, the relative contributions to the spectra of 
neutrinos that interact in the detector are shown. The 
integral for the combined yield from v^, v T , 9^, and v r 
clearly dominates. Further, since the differential cross 
section favors large T p , and since T p ~ E%/M p , the cor- 
responding proton recoil kinetic energy spectrum will be 
much harder, so that they will be even more dominant 
above a realistic detector threshold. 

Since the struck protons in^ + p— > v + p have a rela- 
tively low-energy recoil spectrum, and since realistic de- 
tectors have thresholds, it is crucial to consider the pro- 
ton spectrum in detail, and not just the total yield of 
neutrinos that interact. 



IV. PROTON RECOIL SPECTRUM 

The elastically-scattered protons will have kinetic en- 
ergies of a few MeV. Obviously, these very nonrelativis- 
tic protons will be completely invisible in any Cerenkov 
detector like Super-Kamiokande. However, such small 
energy depositions can be readily detected in scintilla- 
tor detectors such as KamLAND and Borexino. We first 
consider the true proton spectrum, and then in the next 
Section, we consider how this spectrum would appear in 
a realistic detector. 

The true proton spectrum (for one flavor of neutrino) 
is given by 

(T P ) = C [ dE v f(E v ) {E V1 T P ) , (4.1) 

where f{E v ) is a normalized Fermi-Dirac sp ectr um and 
the differential cross section is given by Eq. ( [2.l|) . For a 
given T p , the minimum required neutrino energy is 

T p + ^T p {Tp + 2MJ [M^r 
(E u ) min = ~ y — - — . (4.2) 

The normalization constant C is determined by Eq. ( |3.l| ), 
as the integral of Eq. ( |4.l|) over all T p without the C 
factor is (a). 

In Fig. |3J, we show da/dT p weighted by normalized 
Fermi-Dirac distributions of various temperatures, for a 
single neutrino flavor. Throughout this paper, we refer 
to the u e (T = 3.5 MeV), v e (T = 5 MeV), and the com- 
bined Vp, v T , v^, and v T (T — 8 MeV) flavors. Since we 
know that there are neutrino oscillations, this language 



is somewhat incorrect. However, our results are totally 
insensitive to any oscillations among active neutrinos or 
antineutrinos (since this is a neutral-current cross sec- 
tion), and also to oscillations between active neutrinos 
and antineutrinos (since the cross section is dominated 
by the c 2 A terms). Thus when we refer to the v e flavor, 
we mean "those neutrinos emitted with a temperature 
T = 3.5 MeV, whatever their flavor composition now," 
etc. 

The true proton spectra corresponding to the various 
flavors are shown in Fig. ^. As seen in the hgure, the 
contributions of v e and v e are quite suppressed relative 
to the sum of v^, v r , v^, and v r . 

V. QUENCHING 
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Low-energy protons lose energy very quickly by ion- 
ization. The energy loss rate dE/dx of nonrelativistic 
particles scales roughly as dE/dx ~ —z 2 //3 2 in this en- 
ergy range JIJ], where z is the particle charge and (3 
its velocity. In contrast to the usual —2 MeV/g/cm 2 
for a minimum-ionizing particle, for few-MeV protons, 
dE/dx 100 MeV/g/cm 2 . Thus even a 10 MeV pro- 
ton will be brought to rest in less than about 0.1 cm. In 
contrast, the hadronic interaction length for the proton 
to scatter from a free or bound nucleon is of order 1 cm 
or larger. Thus the hadronic energy losses can be totally 
neglected; see also Fig. 23.1 of Ref. [[L4j. Because of the 
nonlinear response of the detector to proton recoil ener- 
gies, as we are about to describe, it is important that the 
original proton energy is not shared among two or more 
protons, i.e., from elastic hadronic scattering. 

In a scintillator, there is generally an efficient transfer 
between the ionization loss of a charged particle and the 
detectable scintillation light observed by phototubes. For 
example, in KamLAND, there are approximately 200 de- 
tected photoelectrons per MeV deposited for a minimum- 
ionizing particle like an electron |L5| . 

However, for highly ionizing particles like low-energy 
protons, the light output is reduced or "quenched" rela- 
tive to the light output for an electron depositing the 
same amount of energy. The observable light output 
Eequiv (he., equivalent to an electron of energy E equ i V ) is 
given by Birk's Law |l(| : 

dE e (III I r 



dE/dx 



dx 



1 + k B (dE/dx) 



(5.1) 



where fcg is a constant of the scintillation material, and 
dE/dx is the energy deposition rate, now in MeV/cm 
(and defined to be positive). We assume k B — 0.015 
cm/MeV for KamLAND For small dE/dx, the mea- 
sured light output of a proton is equivalent to that from 
an electron of the same energy. But for dE/dx ~ 100 
MeV/cm, the two terms in the denominator are compa- 
rable, and the light output is reduced. At still higher 
dE/dx, then dE equ i V / dx tends to a constant. Birk's Law 
can thus reflects a saturation effect: once dE/dx is large, 



FIG. 5: The quenched energy deposit (equivalent electron 
energy) as a function of the proton kinetic energy. The Kam- 
LAND detector properties are assumed. 



making it larger does not increase the light output. Ef- 
fectively, if all scintillation molecules along the path of 
the particle are already excited, any further energy de- 
position is not converted to visible scintillation light. 

The proton quenching factor was calculated by inte- 
grating Eq. (5.1) with tables |l7| of dE/dx for protons in 
the KamLAND oil-scintillator mixture lEal: 



E, 



,<?„)= ■ 
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dE 



1 + k B (dE/dx) 



(5.2) 



The observed energy in terms of the proton kinetic en- 
ergy is shown in Fig. [5]. Thus the proton quenching factor 
(Eequiv /T p ) is thus roughly 1/2 at 10 MeV, 1/3 at 6 MeV, 
1/4 at 3 MeV, and so on. The detector response is nonlin- 
ear, though in well-understood way. A similar calculation 
using a particles recovered the quenching factor of ap- 
proximately 1/14 noted in Ref. |L5j . Since the energy de- 
position scales roughly as dE/dx ~ z 2 /(3 2 , quenching for 
alpha particles is much worse than for protons of the same 
kinetic energy, since dE/dx is approximately 4 x 4 = 16 
times larger. Our results for the proton quenching fac- 
tor are in good agreement with direct measurements in 
a variety of scintillators jl6], |ll| . 

Using the quenching function shown in Fig. |5|, we can 
transform the true proton spectrum shown in Fig. ^ into 
the expected measured proton spectrum, shown in Fig. ||. 
If the quenching factor were a constant, it would simply 
change the units of the T p axis. However, it is nonlinear, 
and reduces the light output of the lowest recoils the 
most. This increases the effect, shown in previous figures, 
that the measurable contribution from v e and v e is highly 
suppressed relative to the sum of v^, v T , v^, and v T . 

As shown, quenching distorts the spectra according to 
a known nonlinear function. It also reduces the num- 
ber of events above threshold. The anticipated thresh- 
old in KamLAND is 0.2 MeV electron equivalent energy 
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FIG. 6: Analogous to Fig. ^| the struck proton spectrum 
for the different flavors, but with quenching effects taken into 
account. In order of increasing maximum kinetic energy, the 



contributions from u e , u e , and the sum of ty, v T 



and 



v T are shown with dashed lines. The solid line is the sum 
spectrum for all flavors. We assume a 1 kton detector mass 
for KamLAND. 



FIG. 7: The number of events above threshold in KamLAND 



as a function of E t hr in electron equivalent energy E e , 



In 



order of increasing maximum kinetic energy, the contributions 
from v e , y e , and the sum of u^, v T , P M , and v T are shown with 
dashed lines. The solid line is for the sum of all flavors. We 
assume a threshold of 0.2 MeV electron equivalent energy in 
KamLAND. 



(strictly speaking, KamLAND and Borexino have some- 
what higher target thresholds of approximately 0.28 and 
0.25 MeV, set by background rates; over the short du- 
ration of the supernova pulse, much higher background 
rates can be tolerated). With the expected proton 
quenching, this corresponds to a threshold on the true 
proton kinetic energy of 1.2 MeV. The number of events 
above this threshold for each flavor appears in Table | 
Above an electron equivalent threshold of 0.2 MeV, the 
neutrino-proton elastic scattering yields from v e and v e 
are quite small. Thus the measured proton spectrum will 
primarily reflect the shape of the underlying Fermi-Dirac 
spectrum for the sum of v^, v r 9^, and 9 T . Of course, 
this has been convolved with both the differential cross 
section (which gives a range of T p for a given E„), and 
also the effects of quenching. However, as we will show, 
the properties of the initial neutrino spectrum can still be 
reliably deduced. The numbers of events above a given 
electron equivalent threshold are shown in Fig. ^. 



VI. BACKGROUNDS 

In this Section, we consider several backgrounds to the 
signal of neutrino-proton elastic scattering from super- 
nova v^, v T , 9^, and v T . 

The first is neutrino-proton elastic scattering from v e 
and 9 e . As shown above, in particular in Figs. ^ and |t], 
this contribution is minimal above the expected thresh- 
old. We assume that it can be statistically subtracted 
using knowledge of v e and 9 e temperatures measured in 
charged-current reactions and do not consider it further. 



The second comes from a variety of other charged- 
current supernova neutrino signals in the detector. As 
noted above, approximately 310 events are expected from 
v e + p — > e + + n These events can easily be identi- 
fied by the tight coincidence (roughly a few times 10 cm 
in position, and 0.2 ms in time) in the detection of the 
high-energy (about 20 MeV) positrons and the subse- 
quent neutron captures on protons (2.2 MeV gamma). 
There are also charged-current reactions of v e and v e on 
12 C, proceeding almost exclusively to the ground states of 
12 N and 12 B, respectively. About 10 events are expected, 
and possibly a few times more if oscillations effectively 
swap spectra |lq]. However, these events can be iden- 



TABLE I: Numbers of events in KamLAND (1 kton mass 
assumed) above the noted thresholds for a standard super- 
nova at 10 kpc, for the separate flavors or their equivalents 
after oscillations. Oscillations do not change the number of 
neutrinos at a given energy, and the neutral-current yields 
are insensitive to the neutrino flavor. Equipartition among 
the six flavors is assumed (see the text for discussion). The 
thresholds are in electron equivalent energy, and correspond 
to minimum true proton kinetic energies of and 1.2 MeV. 
As discussed in Section II, weak magnetism corrections are 
not included. 
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tified by the subsequent 12 N and 12 B beta decays, with 
lifetimes of order 10 ms and electron endpoints of or- 
der 15 MeV. The total yield from neutrino-electron elas- 
tic scattering (technically, mixed charged- and neutral- 
current) is expected to be about 20 events. We assume 
that these events can be statistically subtracted from the 
spectrum, or that particle identification by pulse-shape- 
discrimination (PSD) will be possible. 

The third comes from other neutral-current supernova 
neutrino signals in the detector. The best-known is the 
superallowed neutral-current excitation of the 15.11 MeV 
state in 12 C, which decays by gamma emission. About 60 
events are expected, and they will be easily identified by 
their narrow spectrum at 15.11 MeV jlq]. There are also 
inelastic neutral-current excitations oP^C that decay by 
particle emission. The yield from all channels that emit 
a proton is about 45 events, using the cross sections and 
branching ratios of Ref. ]nj. Some fraction, probably 
most, will be above threshold and will add to the signal 
of low-energy protons. Though the proton spectra for 
these reactions have not been published, we assume that 
their contribution will be known and can be subtracted. 
The yield from all channels that emit a neutron is about 
20 events, again using the results of Ref. ]l9[ ]. These 
neutrons will be captured, giving 2.2 MeV gamma rays. 
These inelastic neutral-current reactions with proton and 
neutron emission have not previously been recognized as 
supernova neutrino detection channels. Finally, since a 
50 MeV neutrino corresponds to a wavelength of 4 fm, 
about the diameter of a carbon nucleus, there can also 
be coherent neutral-current scattering of the whole nu- 
cleus [^of . The number of events is very large, coinciden- 
tally as large as the total neutrino-proton elastic scatter- 
ing yield (neglecting the detector threshold). However, 
the expected recoil kinetic energies are of course about 12 
times smaller than for free protons. Additionally, since it 
is spin-independent vector scattering, the smallest recoil 
energies are favored. (In contrast, neutrino-proton elas- 
tic scattering is spin-dependent, and the proton spin is 
flipped in the scattering.) The recoil carbon ions will be 
very heavily quenched, and so this signal is unobservable 
in a detector like KamLAND. 

The fourth comes from cosmic-ray induced detector 
backgrounds. Because it is located deep underground, 
the muon rate in KamLAND scintillator is only about 0.3 
Hz, and so all muon-related backgrounds are very small 
over the short duration of the supernova burst |p_ 5|| . 

The fifth and most serious comes from low-energy ra- 
dioactivities in and around the detector. Normally, these 
are not a concern for relatively high-energy supernova 
events. However, here we are considering signals down 
to about 0.2 MeV detected energy, where many different 
radioactive backgrounds contribute. At present, Kam- 
LAND is configured to detect few-MeV reactor antineu- 
trinos via the coincidence between the prompt positrons 
and the delayed neutron captures, and low singles back- 
grounds above 0.2 MeV are not required. Published 
data on the KamLAND background spectrum are not 



yet available. However, if KamLAND is to be eventu- 
ally used for detecting Be solar neutrinos by neutrino- 
electron scattering, then the background in this energy 
range will have to be reduced to about 10 -3 Hz, the rate 
of solar neutrino events expected (similar considerations 
hold for Borexino). For the supernova signal discussed in 
this paper, a much larger background rate of about 1 Hz 
could be tolerated. This rate is set by the consideration 
of being much less than (300 events/10 s) = 30 Hz. 

Therefore, in what follows we consider just the main 
signal, and neglect backgrounds. 

VII. PROTON SPECTRUM FITS 

In this Section, we show how the measured proton 
spectrum can be used to separately determine the to- 
tal energy of the v^, is T , v^, and v r neutrinos and their 
time-averaged temperature. The total number of de- 
tected events is proportional to the portion of the to- 
tal binding energy carried away by these four flavors, 
and we denote this by E tot (note that this is not the 
total binding energy Eb)- For a standard supernova, 
E tot = 4(£ B /6) = 2/3E B ~ 2 x 10 53 ergs. We denote 
the temperature of these four flavors by T . If only the 
total yield were measured, as for most neutral-current 
reactions, there would be an unresolved degeneracy be- 
tween E tot and T, since 

N~E tot ^-. (7.1) 

Note that for er ~ E™, then (a) ~ T n . For v + d 
v + p + n in SNO, for example, a ~ E 2 , so N ~ E tot T. 
Thus for a given measured number of events, one would 
only be able to define a hyperbola in the plane of E tot and 
T. The scaling is less simple here because of threshold 
effects, but the idea is the same. 

Here we have crucial information on the shape of the 
neutrino spectrum, revealed through the proton spec- 
trum. To remind the reader, in most neutral-current re- 
actions there is no information on the neutrino energy, 
e.g., one only counts the numbers of thermalized neutron 
captures, or measures nuclear gamma rays (the energies 
of which depend only on nuclear level splittings). 

In this Section, we perform quantitative tests of how 
well the parameters E tot and T can be determined from 
the measured proton spectrum. (We did also investigate 
the effects of a chemical potential in the Fermi-Dirac dis- 
tribution, but found that it had little effect. This is sim- 
ply because the cross section is not rising quickly en oug h 
to see the tail of the thermal distribution in detail plf ■) 
Of course, if the distance to the supernova is not known, 
then we are effectively fitting for E tot / D 2 . 

We performed Monte Carlo simulations of the super- 
nova signal in KamLAND and made chi-squared fits to 
determine E tot and T for each fake supernova. To per- 
form the fits, we started with an "ideal" spectrum, as 
described by the integral: 




2 



FIG. 8: Example spectra with different values of E tot and 
T, all chosen to give the same number of events above an 
electron equivalent threshold of 0.2 MeV (true proton energy 
1.2 MeV) in KamLAND. Though not shown in this figure, the 
spectrum above 2 MeV is included in our analysis. At the 0.2 
MeV point, from left to right these correspond to (E tot ,T) 
= (4.2, 6), (2.0, 8), (1.4, 10), respectively, with E tot in 10 53 
ergs and T in MeV. (Thus with the standard E to t = 2 x 10 53 
ergs, the number of events above threshold with T = 6 MeV 
is 2.0/4.2 times the number with T = 8 MeV; with T = 10 
MeV, it is 2.0/1.4 times the number with T = 8 MeV). 

fdN\ f°° 
"PF = C dT;G(T;;T p ,6T p ) 

\ al p/ ideal Jo 

POD J 

x / dE v f{E„)—{E u ,T' p ),{1.2) 

J(E v ) mln dT p 

where the inner integral is as in Eq. fl4.1p , with the ad- 
dition that quenching corrections are applied to T' v after 
convolution with f{E v ). For the Gaussian energy resolu- 
tion G(Tp] T p , ST P ), we used ST P = 0. W(l MeV) g|. 
Because the proton spectrum has already been smeared 
by the neutrino spectrum and the differential cross sec- 
tion, the Gaussian energy resolution has only a minor 
effect. The norm alization constant C is given by compar- 
ison to Eq. ( |3.l[ ). Example spectra are shown in Fig. |^. 

Using (dN/ dTp) id ca i, we binned the spectrum by the 
following integral: 

Ni= dT p[^r) ( 7 - 3 ) 

J(.E mtn )i \ al pJ ideal 

where Ni is the number of events in bin i, and (E m i n )i 
and (E max )i are the minimum and maximum energies for 
bin i. Eight bins of variable width were used, chosen to 
contain roughly the same number of expected events per 
bin. For a chosen E tot and T, this was the starting point 
of our Monte Carlo (and the bin boundaries were kept 
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FIG. 9: Scatterplot of 10 3 fitted values, in the E tot and T 
plane, for the labeled "true" values, where E tot is the total 
portion of the binding energy carried away by the sum of 
Vfi, v T , Vfi, and u T , and T is their temperature. The values 
of E tot and T were chosen such that the numbers of events 
above threshold were the same. The measured shape of the 
proton spectrum breaks the degeneracy between these two 
parameters. Without that spectral information, one could 
not distinguish between combinations of E tot and T along the 
band in this plane that our three example regions lie along. 

fixed). For each fake supernova, we sampled the number 
of events in each of these bins according to the appropri- 
ate Poisson distributions. The resulting spectrum was as 
one might obtain from a single supernova, given the finite 
numbe r of events expected. We then varied E tot and T 
in Eq. ( |7.3| ) until the values that best fit the fake spectral 
data were determined. For a given set of assumed E tot 
and T, this procedure was repeated many times. The 
distributions of the final E tot and T thus reveal the ex- 
pected errors on fitting E tot and T for a single real future 
supernova. 

Three examples are shown in Fig. ^, where one can see 
that E tot and T can each be determined with roughly 
10% error. These errors scale as 1/yN, where N is the 
total number of events (i.e., if one imagines a detector of a 
different mass or a different assumed supernova distance) . 
If the distance were completely uncertain, one would not 
be able to determine E tot . However, after marginalizing 
over the unknown E tot (i.e., projecting these scatterplots 
onto the T axis), one would still obtain a good measure- 
ment of T. 



VIII. DISCUSSION AND CONCLUSIONS 

We have shown that neutrino-proton elastic scattering, 
previously unrecognized as a useful detection reaction for 
low-energy neutrinos, in fact has a yield for a supernova 
comparable to v e + p — > e + + n, even after taking into 
account the quenching of the proton scintillation light 
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and assuming a realistic detector threshold. 

In addition, the measured proton spectrum is related 
to the incident neutrino spectrum. We have shown ex- 
plicitly that one can separately measure the total energy 
and temperature of v^, v T , D^, and v T , each with uncer- 
tainty of order 10% in KamLAND. This greatly enhances 
the importance of detectors like KamLAND and Borex- 
ino for detecting supernova neutrinos. 

For Borexino, the useful volume for supernova neutri- 
nos is 0.3 kton, and the hydrogen to carbon ratio in the 
pure pseudocumene (C9H12) is 1.3 : 1 [ p2[ , so there are 
about 4.7 times fewer free proton targets than assumed 
for KamLAND. However, the quenching is less in pure 
scintillator (KamLAND is about 20% pseudocumene and 
80% paraffin oil and the errors on E tot and T scale 

as l/\/~N, so that the precision in Borexino should be 
about 20% or better. 

Other techniques for bolometric measurements of su- 
pernova neutrino fluxes have been studied. Detectors 
for elastic neutral-current neutrino scattering on elec- 
trons J2|| and coherently on whole nuclei [^0| have been 
discussed, but never built. If neutrino oscillations are ef- 
fective in swapping spectra, then the temperature of the 
"hot" flavors may be revealed in the measured positron 
spectrum from D e + p — > e + + n; two recent studies have 
shown very good precision (< 5%) for measuring the tem- 
peratures and the total binding energy (24|, ^H). However, 
they assumed exact energy equipartition among the six 
neutrino flavors, whereas the uncertainty on equiparti- 
tion is at least 50% [[ll|. Nevertheless, under less re- 
strictive assumptions, this technique may play a comple- 
mentary role. Finally, since for different cross sections, 
the neutral-current yields depend differently on temper- 
ature, comparison of the yields may provide some infor- 
mation [p6 |. However, there are caveats. In neutrino- 
electron scattering, the neutrino energy is not measured 
because the neutrino-electron angle is much less than the 
angular resolution due to multiple scattering. The scat- 
tered electrons, even those in a forward cone, sit on a 
much larger background of v e + p — > e + + n events, so 
it is difficult to measure their spectrum |27 ; also, their 
total yield is only weakly dependent on temperature. At 
the other extreme (see Fig. 3 of Ref. ^6|), the yield of 
neutral-current events (2j| on 16 O depends strongly on a 
possible chemical potential term in the thermal distribu- 
tion. 

It is important to note that the detection of recoil 
protons from neutron-proton elastic scattering at several 
MeV has been routinely accomplished in scintillator de- 
tectors (see, e.g., Ref. @). Since both particles are mas- 
sive, the proton will typically take half of the neutron en- 
ergy. This reaction provides protons in the same energy 
range as those struck in neutrino-proton elastic scatter- 
ing with E v ~ 30 MeV. This is a very important proof 
of concept for all aspects of the detection of low-energy 
protons. 

Though low-energy backgrounds will be challenging, it 
is also important to note that the background require- 



ments for detecting the supernova signal are approxi- 
mately 3 orders of magnitude less stringent than those 
required for detecting solar neutrinos in the same en- 
ergy range (taking quenching into account for our signal) . 
Borexino has been designed to detect very low-energy so- 
lar neutrinos, and KamLAND hopes to do so in a later 
phase of the experiment. 

These measurements would be considered in combi- 
nation with similar measurements for v e and v e from 
charged-current reactions in other detectors. Separate 
measurements of the total energy and temperature for 
each flavor will be invaluable for comparing to numerical 
supernova models. They will also be required to make 
model-independent studies of the effects of neutrino os- 
cillations. If the total energy release Eb in all flavors has 
been measured, then 



E 



B — T 



3 GM% S 
5 Rns 



(8.1) 



thus allowing a direct and unique measurement of the 
newly-formed neutron star properties, principally the 
mass Mns p9|]. 
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